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Neutrino and axion production in neutron stars occurs 
mainly as bremsstrahlung from nucleon-nucleon (NN) scat- 
tering. The energy radiated via neutrinos or axions is typi- 
cally very small compared to other scales in the two-nucleon 
system. The rate of emission of such "soft" radiation is di- 
rectly related to the on-shell NN amplitude, and thereby 
to the NN experimental data. This facilitates the model- 
independent calculation of the neutrino and axion radiation 
rates which is presented here. We find that the resultant rates 
are roughly a factor of four below earlier estimates based on 
a one-pion-exchange NN amplitude. 



Neutron stars are believed to be born during a su- 
pernova explosion with an interior temperature T of or- 
der 60 MeV. The subsequent evolution of the hot and 
dense compact star is characterized by a rapid early cool- 
ing phase followed by a, significantly slower, late-time 
cooling phase |Q. During both of these phases neutri- 
nos are an important source of energy loss. Thermal 
evolution of neutron stars is largely driven by neutrino 
bremsstrahlung reactions such as 



NN -► NNvv 



ripe v e . 



(1) 



In the first part of the nascent neutron star's life it 
evolves by diffusion of trapped neutrinos. This cools the 
interior to T ~ 1 MeV on a time scale of a few seconds . 
The resultant intense neutrino emission is thought to play 
a central role in both the supernova mechanism and r- 
process nucleosynthesis. The high neutrino luminosity 
is fueled by the reactions (|l|), which compete with the 
annihilation e + + e~ — * vv in degenerate matter. At 
later times, the neutron star enters a period of slower 
thermal evolution, during which the emitted neutrinos 
free stream. This occurs because the neutrino mean-free 
path becomes long when T < 1 MeV. The time scale for 
this long-term cooling of the dense, degenerate, neutron- 
rich, inner core thus depends crucially on the neutrino 
emissivity, which is, again, dominated by the reactions 
(|l|). Observational constraints on this late-time portion 
of the neutron star's evolution will improve as X-ray ob- 
servatories such as Einstein, EXOSAT and ROSAT @] 
gather pulsar data which gives information on the surface 
temperatures of these stars. The challenge for theorists 
is to improve the models of both the early- and late-time 
cooling of the neutron star. 



The emissivities which are key ingredients in these sim- 
ulations are dominated by the reactions ([!]) . Despite their 
central role in neutron-star dynamics these reactions have 
received relatively little attention. Pioneering work was 
done by Friman and Maxwell ||] who computed the reac- 
tion rates using a nucleon-nucleon amplitude due only to 
a single pion exchange (henceforth we refer to such calcu- 
lations as the "OPE approximation"). Recently, many- 
body effects, in particular the suppression due to multiple 
scattering (the Landau-Pomeranchuk-Migdal, or LPM, 
effect ^|), have been shown to be important at temper- 
atures T> 5 MeV [|-§. However, the work of Ref. [§ is 
still the state of the art treatment of the NN interactions 
which occur during the reactions (|l|) (see also Ref. ||). 
In this letter we follow an approach reminiscent of that 
used in soft-photon calculations jio) , and relate the rate 
of production of soft-neutrino radiation in NN scattering 
to the on-shell NN scattering amplitude. This yields a 
calculation we present as a "benchmark" , which accounts 
for the two-nucleon dynamics in a model-independent 
way. We make no attempt to account for many-body 
effects, although they are undoubtedly important in the 
star. We identify the density and temperature range over 
which our results are valid and show that, although lim- 
ited, it is of interest to both supernova and neutron star 
physics. In contrast, the full evaluation of the rates for 
the reactions ([!]) in a strongly-coupled medium is a com- 
plicated problem. Therefore, of necessity, most solutions 
to it will be model-dependent. Thus, we see our results 
as providing a model-independent foundation on which 
future work that assesses the role of many-body effects 
can build. 

Our main focus in this work will be the neutrino emis- 
sivity from NN — > NNvv. However, axions (if they 
exist) couple, like neutrinos, to the nucleon spin. There- 
fore, in computing the neutrino emissivity one obtains 
the axion emissivity from NN — > NNa as a welcome by- 
product ||. This is useful because one important con- 
straint on the axion mass comes from considering the role 
of axionstrahlung in the dynamics of SN1987A In- 
deed, if the rate for the reaction NN — * NNa is too high 
then the supernova dynamics is completely changed, and 
the successful "standard" picture of the supernova is de- 
stroyed. Thus, one can constrain the axion coupling, and 
hence the axion mass, by demanding that axion radiation 
did not make too large a contribution to the energy loss 
from SN1987A. 

Neutrino and axion emissivities: We begin by explic- 
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itly calculating the emissivity due to NN — > NNvv. The 
vv coupling to non-relativistic baryons at low energies is 
given by the Lagrange density 

G 



(2) 



where f = vj^^l — r y 5 )v is the leptonic current, Gf — 
1.166 x 10~ 5 GeV~ 2 , N is the nucleon field, and c v and c a 
are the nucleon neutral-current vector and axial-vector 
coupling constants. Some Feynman diagrams for the 
bremsstrahlung process are shown in Fig. |l|. 
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FIG. 1. Feynman diagrams for the bremsstrahlung process. 
The radiation is represented by the dashed line, and nucleons 

(2) 

by solid lines. Tatjv is the AW transition matrix and J^ 5 is 
a two-body axial current. 

The incoming (outgoing) nucleon momenta are labeled 
PijP2 (P3,p 4 )- The dashed line represents radiation — 
a neutrino-anti-neutrino pair in this case — which carries 
energy lo and momentum q. In general we are interested 
in cases where the radiated energy is small compared 
to the incoming nucleon energy. In the limit lo —> the 
amplitudes corresponding to diagrams (a) and (b) in Fig. 
1 are dominant, as they contain pieces proportional to 
1/lu. On the other hand, the contributions from the re- 
scattering diagram Fig. 1(c), and from meson-exchange 
currents such as Fig. 1(d), remain finite in the lo — > 
limit. Thus, for the reaction nn — > nnvv the matrix 
element can be written as 
Gf 



M = 2 



- l^(p'\[T NN ,r,}\p) + 0(^), (3) 



where p (p') is the initial (final) relative momentum of 
the two-nucleon system. We refer to results which retain 
only this leading term, of 0(lo~ 1 ), in A4 as "true in the 
soft-neutrino approximation (SNA)" . In general the NN 
T-matrix appearing in Eq. (y), Tnn, will be half off- 
shell []. But, in the SNA we can take Tnn to be the on- 



shell NN amplitude. We can also neglect the difference 
between the magnitude of the initial and final-state rela- 
tive momenta. We expect these approximations to break 
down when lo ~ m OT , since sets the scale for varia- 
tions of Tn n in the off-shell direction ^ So, in the SNA, 
the NN interaction is described by the on-shell T-matrix 
Tatjv, evaluated at a center-of-mass energy which, for 
reasons of symmetry, is chosen to be (p 2 + p' 2 ) / (2M) (M 
is the nucleon mass). This T-matrix can be constructed 
from phase shifts deduced from NN scattering data |L2j . 
Note that the OPE approximation used in most previous 
calculations involves substituting Vope, the one-pion- 
exchange potential, for Tnn in Eq. (0). Meanwhile, 
is the vertex which couples the radiation to the nucleons. 
For vv radiation T M follows straight from Eq. (0). Only 
its three-vector part contributes to M. at 0{Lo^ r ). Equa- 
tion (||) then gives us a model-independent result for M , 
which is correct in the SNA. 

If only two-body collisions are taken into account then 
the neutrino emissivity from a neutron gas is given by 
Fermi's golden rule 



f - 



(2^)32^! (2tt) 3 2w ; 



-(2ir) 4 6(E m -E fn ) 



LO 5 3 (pi, 



Pfr. 



d 3 pi 



i=1..4 v ; 



T X - £|A4| 2 , (4) 



where T = /x/ 2 (l - / 3 )(1 - / 4 ), with f % = 1/(1 + 
exp [Ei — fj,i)/T) being the Fermi-Dirac distribution 
function for the nucleons, and s = 4 the symmetry fac- 
tor accounting for identical nucleons. The spin-summed 
square of the matrix element can be factored into leptonic 
and hadronic tensors, and then represented by 



Tr (lH j ) Hi 



(5) 



The trace over the lepton tensor is easily evaluated. Fur- 
ther, since we are interested in soft radiation, we may 
safely ignore q in the momentum delta function |jj . This 
allows us to directly integrate the leptonic trace over neu- 
trino angles to obtain 



J dCli J dfl 2 TR (klj) = 



8 (4tt) 2 loilo 2 S. 



t-,3 



(6) 



Therefore, only the trace of the hadronic tensor Tiij con- 
tributes to the emissivity, and so we define a scalar func- 
tion, 



As used here, it should involve a sum over the allowed 
partial-waves of the AW system. This, together with the 
factor of two in front of the matrix element, accounts for the 
exchange graphs which must be included in M. 



2 At very low relative momenta the scale of breakdown is set 
by the AW scattering length, since that gives the variation in 
the on-shell direction. However, ajvjv does not really play a 
role here, since typical nucleon momenta in neutron stars are 
at least 100 MeV. 
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(7) 



which is called the dynamical spin structure function of 
the medium. It is related to the vv emissivity via: 



r 2 r 2 1 

(j F C a 1 
167T 4 30 



doj uj 6 S a {to) , 



(8) 



where u> is the total energy of the emitted vv pair. 

In the two-body approximation considered here we 
evaluate Tin using Eqs. (|) and (|). For the case of 
emission from the nn system, only the spin-triplet two- 
nucleon state contributes, and the trace is: 

n u = l6-^ \(lM' s ,p'\[S t ,T NN }\p,lM s )\ 2 , (9) 



M a M' 3 

where Si is the total spin of the two-nucleon system. It is 
straightforward to generalize this formula for Tin to the 
np case, although the NN spin singlet then contributes. 
(The np case, and the failure of the OPE approximation 
there, is discussed in Ref. (p|.) From Eqs. (||) and (Q) 
we can calculate S&, and thus the vv emissivity. 

The emission of any radiation which couples to the 
nucleon spin will be described by the same function S&. 
Thus, as mentioned above, with S a in hand we may de- 
rive the axion emissivity £ a . The effective theory for 
axion-nucleon interactions is described by the Lagrange 
density £ allI1 = — g ann aN~f 5 N, where a is the axion field, 
and g ann = 10~ 8 (rn a /l cV) is the effective axion coupling 
(m a is the axion mass) |14 . The calculation of the ax- 



ion emissivity in this effective theory is analogous to the 
above calculation of the neutrino emissivity, and yields 



ffann 1 
167T 2 Af 2 3 



dui u> S a (uj) 



(10) 



Before proceeding to our results we note that S& can 
be defined in a much more general way, where it describes 
the response of a many-body system to an external spin- 
dependent perturbation. Equations (j|) and ([W) remain 
true if this definition is adopted. To obtain Eq. ( 71) for S a 
in this general case one takes the long-wavelength limit 
of the leading term in the density expansion. 

Results & Discussion: We present results for the dy- 
namic spin structure function S a (uj), since it includes 
the density, temperature, and nuclear dynamics depen- 
dence of the neutrino and axion emissivities. During the 
evolution of neutron stars, one encounters varying de- 
grees of nucleon degeneracy, with (X n /T ~ 1 at birth, but 
Hn/T ^> 10 at late times. Earlier investigations have 
shown that analytic approximations to the phase-space 
integrals in Eq. (Q) work poorly at intermediate degen- 
eracy uW. Therefore, in this work these integrals are 



all performed numerically. In order to investigate the 
effect of our model-independent treatment of the NN in- 
teraction we plot the ratio R a (u) = S% NA (u)/S r a ci (w), 
where S^ na (uj) is calculated as described above. The 
denominator, S 1 ^ 1 (ui), is the dynamic spin structure func- 
tion found when a hadronic tensor trace of the form 
Tin = c/uj 2 is inserted into Eq. (Q). We adjust the con- 
stant c so that when S l a ot is employed in Eq. (||) the neu- 
trino emissivity thereby obtained is equal to that found if 
the full OPE matrix element is used in evaluating Tin []. 
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FIG. 2. A plot showing the ratio R a = S^ NA /S r f for nn 
pairs in neutron matter at a density of 0.16 fm -3 . The region 
where multiple scattering suppression due to the LPM effect 
is negligible (u> > 7, see below), and the SNA is expected to 
be valid (u < m^) is delineated by thicker curves. The regions 
probed by the neutrino and axion emissivities are different, 
as indicated on the plot. 

Figure ^ shows the resulting ratio R<j(lo) for neutron 
matter at a range of temperatures and a baryon density 
equal to the nuclear saturation density. The results are 
plotted as a function of the dimensionless ratio oj /T (note 
S a (oj) has significant strength only for lo/T < 15). 

The most striking feature of the results is that the one- 
pion-exchange approximation significantly overestimates 
the rate for neutrino (or axion) production. The large re- 
duction in the response functions over those obtained in 
the OPE approximation occurs for two reasons. Firstly, 
one-pion exchange over-estimates the strength of the NN 
tensor force, and so even replacing the Vope employed 
previously with the full Vnn would lead to a reduction 
in S a . Secondly, the unitarity of our NN T-matrix leads 
to an NN amplitude which is, in general, significantly 
smaller than that found in the OPE approximation, and 



3 We could have followed Refs. |]|ll| and adopted a reference 
Scf(u) in which the full OPE-approximation matrix element 
is replaced by its value in the — > limit. However, this 
is a poor approximation to the actual result for one-pion ex- 
change, since it over-estimates the OPE-approximation emis- 
sivities by as much as a factor of two (see also Refs. |p[|ll]). 
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hence to a much- reduced S a . That the OPE approxima- 
tion does such a poor job in describing the NN dynamics 
should come as no surprise. At the energies of interest 
here the NN interaction is intrinsically non-perturbative, 
and so replacing the full AW T- matrix by Vope gives 
only a crude estimate of these rates. 

Fig. ^ may be used to infer where the SNA breaks 
down. Recall that S vei is constructed using a hadronic 
tensor Hij proportional to 1/w 2 . In fact this part of 
the hadronic response is the only piece of the Hij cal- 
culated in the SNA that can be trusted. Therefore, any 
^-dependence of R a (oj) represents an effect in S% NA (w) 
coming from physics beyond the SNA. Viewing Fig. |2| in 
this light, and considering the constraint u> < m^, implies 
that at p = po the SNA works well for T < 10 MeV. 

The suppression of the axial response seen in Fig. || 
translates into a corresponding diminution of axion and 
neutrino emissivities. Let us define the ratio Rg u = 
£^ A /£°f E . Here, £°f E is the emissivity found in the 
OPE approximation, as calculated in Ref. 0. (A similar 
ratio of axion emissivities is approximately equal to Rs u 
in the domain of validity of the SNA.) The ratio Re v is 
displayed in the table below for a range of dcnsitites and 
at temperatures of 1 and 10 MeV. In fact, the tempera- 
ture and density dependence of £^ A is predominantly 
determined by that of the nucleon equilibrium distribu- 
tion functions appearing in Eq. (||), but the ratio Rg v 
changes significantly over the densities and temperatures 
considered because £^f E has much more variation with 
tib and T. We see that, for the range of conditions con- 
sidered here, the SNA gives a rate of emission of soft axial 
radiation which is roughly a factor of four smaller than 
that given by the OPE approximation. 



n B (fm s ) 


R£ (1 MeV) 


(10 MeV) 


0.08 


0.29 


0.27 


0.16 


0.24 


0.23 


0.48 


0.16 


0.16 



Disclaimers: As mentioned previously, our calcula- 
tion makes no attempt to include many-body effects. 
These will certainly be important in some regimes of tem- 
perature and density. For instance, it is claimed that 
in-medium modifications of the pion, attributed to the 
many-body nature of the problem, strongly affect the 
emission rates Such effects are outside the scope of 
this work. However, as stated earlier, we expect the LPM 
effect to strongly reduce the response function J|-^]. In 
particular, it will suppress the emission of radiation with 
lo < 7, where 7 is the nucleon quasi-particle width at 
the Fermi surface. This LPM-effect limit on the validity 
of our calculation is indicated in Fig. with the value 
of 7 taken from Ref. jjj . Figure || suggests that LPM- 
suppression will affect the emissivity if T > 10 MeV. Note 
that the LPM effect and the use of the SNA both signifi- 
cantly reduce the rate of emission of soft axial radiation. 



Several microphysical ingredients play a role in the 
thermal evolution of neutron stars, and so it is difficult 
to state precisely how the results obtained in this work 
will affect observable aspects of neutron star evolution. 
However, our results imply that NN bremsstrahlung is 
less important during the star's infancy than previously 
thought. In addition, the reduction of the axion emissiv- 
ity we have found will, presumably, weaken the axion- 
mass bound obtained from SN1987A. However, these 
comments are subject to the caveat that through large 
regions of the infant neutron star the temperature is high 
enough to invalidate the SNA. 

On the other hand, for late-time cooling the tempera- 
ture is small, and our results are applicable everywhere in 
the star. However, in this regime the modified URCA re- 
action, nn — > npe~v ei is significantly more efficient than 
the pair process considered here || . In degenerate matter 
this charged-current reaction does not produce soft radi- 
ation, since the typical change in the energy of the NN 
system is of order the electron chemical potential, i.e. 
about 100 MeV. Nevertheless, the results presented here 
suggest that large corrections to the modified URCA rate 
calculated in the OPE approximation will occur when a 
better model of the NN amplitude is used. 

Conclusion: Finally, we reiterate that none of these 
disclaimers modify the two central conclusions of this pa- 
per: that the soft-neutrino approximation gives a model- 
independent result for the emissivity due to the reactions 
NN -> NNvD and NN -> NNa; and that these emis- 
sivities are much smaller than those found when one-pion 
exchange is used as the NN amplitude. 
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